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Summary: The field of exoplanetary science is one of the most rapidly growing areas of 
astrophysical research. As more planets are discovered around other stars, new techniques 
have been developed that have allowed astronomers to begin to characterise them. Two of 
the most important factors in understanding the evolution of these planets, and potentially 
determining whether they are habitable, are the behaviour of the winds of the host star and 
the way in which they interact with the planet. The purpose of this project is to reconstruct 
the magnetic fields of planet hosting stars from spectropolarimetric observations, and to use 
these magnetic field maps to inform simulations of the stellar winds in those systems using 
the Block Adaptive Tree Solar-wind Roe Upwind Scheme (bats-r-US) code. The bats-r- 
US code was originally written to investigate the behaviour of the Solar wind, and so has 
been altered to be used in the context of other stellar systems. These simulations will give 
information about the velocity, pressure and density of the wind outward from the host star. 
They will also allow us to determine what influence the winds will have on the space weather 
environment of the planet. This paper presents the preliminary results of these simulations 
for the star r Bootis, using a newly reconstructed magnetic field map based on previously 
published observations. These simulations show interesting structures in the wind velocity 
around the star, consistent with the complex topology of its magnetic field. 
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Introduction 

Exoplanetary science is one of the most dynamic and rapidly expanding fields in modern 
astrophysical research. Since the first discovery of a planet orbiting a sun-like star, less than 
2 decades ago [Ql, astronomers have found over fifteen-hundred planets orbiting other star^ 
The early discoveries were of planets of comparable mass to, or more massive than, Jupiter, 
and orbiting at distances much closer to their star than the Earth is from the Sun [[31 |4l |5||. 
Advances in technology, and the increasing length in the temporal baseline of observations, 
have resulted in the discovery of systems that have large planets orbiting at similar distances 
to the gas giants of our Solar system [IH |7l [H. Also, as a result of the exquisite precision 
with which the transits of planets across their host stars can now be measured, the size of 
the smallest discoverable planets has been decreasing, such that the first “super-Earths” have 
now been found [|9l dOl fTTI . 

*As of 6th November, 2014, the current tally stands at 1516 confirmed planets, with a further 3359 candidate planets 
awaiting confirmation, according to the Exoplanet Orbit Database at http://exoplanets.org m. 


As a result of the dramatie progress made in the diseovery of exoplanets, the foeus of 
exoplanetary seienee has shifted to inelude trying to eharaeterise the planets that have been 
found. It will not be long before the first Earth-analogue planets will be diseovered, and the 
hunt for life on those planets will begin. The eharaeterisation of these Earth-like planets will 
be pivotal in deeiding whieh ones to examine more elosely. Determining the habitability of a 
planet based on a large range of variables is vital to ehoosing the most promising target for 
the seareh for life outside our Solar system ffT^ . 

A vital eomponent in determining the habitability of a planet is the infiuenee of the 
planet’s host star. In partieular, the aetivity and winds resulting from the magnetie field of 
an exoplanet’s host-star will play a eritieal role in ensuring that the planet will be able to 
sustain a suffieient atmosphere to ensure that life ean develop and thrive. In our own Solar 
system, Mars stands as an example of the destruetive power winds ean have on a planet’s 
atmosphere. Over the four and a half billion year history of our Solar system, Mars’s onee 
thiek atmosphere has been stripped by the wind of our Sun lfT3l . Earth has been able to 
maintain the atmosphere required to support life due to its proteetive magnetie field. 

It is elear, then, that an understanding of the interaetion between the stellar winds of planet 
host stars and the planets that orbit them will play an important role in the seareh for life 
outside our Solar system. However, this area of researeh is still in its infaney, due to the 
diffieulties involved in either deteeting and measuring stellar winds drreetly or in modelling 
them based on the information we ean obtain about their magnetie field behaviour. 

The first deteetions of Solar-like winds around other stars were made by examining the 
blue-shifted absorption of Eyman-a photons by heated atomie hydrogen gas in the Interstellar 
Medium (ISM) llTdll . This is an indireet method of deteeting winds, and attempts have been 
made to more direetly deteet winds through looking for x-ray emission from the interaetion of 
the winds and the ISM IfTSlI . or radio emission direet from these winds [|T^ . Unfortunately, 
these methods have proven to be two to three orders of magnitude less sensitive than the 
Eya absorption teehnique [fT4ll . As a result of the laek of direet observation, the behaviour of 
stellar winds remains poorly understood. We are thus reliant on the modelling of these stars, 
eonstrained by the observations we ean make, to investigate the behaviour of stellar winds. 

In order to model the winds of other stars, we require information on the behaviour of their 
magnetie fields, and this ean be obtained using Zeeman Doppler Imaging (ZDI) ifTTll . Making 
speetropolarimetric0 observations at different phases of a star’s rotation, ZDI ean reeonstruet 
the topology of the radial, azimuthal and meridional parts of its magnetie field (e.g. IIT9ll20ll ). 

The purpose of this projeet is to take the reeonstrueted radial magnetie fields of planet 
hosting stars from speetropolarimetrie observations, and use these magnetie field maps to 
inform a simulation of their winds using the BATS-R-US eode fl^ . In this paper, we will 
outline our ehosen data for this projeet and target star for initial testing, and briefly diseuss 
the eonstruetion of the magnetie field maps that will be used as input for the stellar wind 
modelling. We will then present some of the preliminary results of the projeet, and outline 
the work to be eompleted in the future. 

^Spectropolarimetry is an observational technique that involves passing the light from a star through a polariser before 
sending it into a high-resolution spectrograph. Since magnetic fields will polarise light from atomic spectra, observing 
polarised spectra then detects the magnetic field of a star. A detailed explanation of this can be found in fTsl . 


Data Sample 


This project will make use of spectropolarimetric data available through our membership 
of the BCoo0 project Il22l . with a focus on investigating solar type, planet hosting stars. The 
first star to be examined for preliminary testing is r Bootis (r Boo), a solar type star with a 
~ 5 Jupiter-mass planet that orbits at a distance of 0.05 AU and with a period of 3.31 days 
01 . Table gives a summary of the basic properties of this star taken from the Exoplanet 
Orbit Database!^ [|3. 

r Boo is an active star, whose magnetic field and wind behaviour has been previously 
studied (e.g. [l23l|24l|25l|2^), and so is an ideal test target. It displays a cycle in its magnetic 
field similar to that of our Sun’s ~ 22 year magnetic cycle, but over a much shorter timescale 
of around 2 years [|25l. Given the number of magnetic observations taken over many epochs, 
this makes it an ideal candidate for investigating the temporal evolution of the wind during 
its magnetic cycle, and the impact this may have on its planet. 

Table 1: Summary of the basic properties of the star r Boo. 


Mass (Solar masses) 

1.34 

Radius (Solar radii) 

1.46 

Stellar usin/ (km s“^) 

14.98 IM 

Stellar Inclination Angle (deg.) 

40 mi 

Apparent Magnitude 

4.5 m 

Effective Temperature (K) 

6387 dm 

Spectral Type 

E6IV mi 


Stellar Magnetic Field Maps 

The spectropolarimetric observations obtained of our target stars are used to map stellar 
magnetic fields using ZDI. The process for reducing the data is based on the latest version 
of the “Echelle Spectra Reduction: an Interactive Tool” software [l32l . called libre-esprit, 
developed by Jean-Erancois Donati of the Observatoire Midi-Pyrenees. This software in¬ 
corporates the optimal extraction of spectra from the CCD images, and uses least squares 
deconvolution (ESD) to co-add spectral lines, significantly increasing the sensitivity of mag¬ 
netic detections in the spectra. This is needed, as the magnetic field signal is only ~ 0.1% 
of the total light from a star ll32l . Since the magnetic field will effect the all spectral lines in 
a similar way, the thousands of spectral lines in each Echelle spectrum can be co-added in 
this way to improve the signal-to-noise of the magnetic field detection. 

The method of ZDI then uses these ESD profiles from observations made at different 
phases in the star’s rotation to reconstruct the map of the magnetic field for each of the 
radial, azimuthal and meridional field component^ Examples of such maps for young, solar 

^The BCool project (http://bcool.ast.obs-mip.fr/l is an international collaboration that is carrying out ongoing observations 
to categorise and understand tbe magnetic fields of cool stars. 

^ http://exoplanets.org/ Data were obtained on the 17th October, 2014. 

^More details on this can be found in ca. 




type stars can be found in ^33[ [34l 1^ . The stellar winds have been shown to be mainly 
governed by the radial component of the magnetic field and so the radial field map 
alone is used as input to the simulations in our work. 

Figure shows the radial magnetic field configuration at the surface of r Boo from 
observations made by Rim Fares and collaborators [l25l in July 2008 using the Naval echelle 
spectropolarimeter on the Telescope Bernard Lyot (TBL) in France. This magnetic field map 
differs slightly from the those published of the same data set previously due to advances in 
the ZDI method. The map shows multiple regions of positive and negative magnetic field, 
with several large positive field spots and a large negative spot around the Northern pole. 
The field is less detailed at the Southern pole due to a lack of information resulting from the 
inclination angle of the star with respect to us. The black line indicates where the field is 
zero. 
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Fig. 1: A plot of the surface magnetic field of r Boo from observation taken in July 2008 
by Fares et al. [l25ll . The black line indicates where the magnetic field is zero. The most 
interesting feature of this plot are the multiple regions of positive and negative magnetic 
field, with a large negative spot around the northern pole, and the multiple, large positive 
regions in the mid-Northern latitudes. No information is available for the far southern region 
of the star as that is pointed away from the Earth. 


Stellar Wind Modelling With bats-r-us 

The wind modelling for our target stars is done using the BATS-R-US code, originally 
developed to simulate the Solar wind [|2T]| . bats-r-us solves the three-dimensional, ideal 
magnetohydrodynamic equations: 
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which describe the conservation of mass, momentum, magnetic flux and energy, respee- 
tively. In these equations, p is the mass density, u is the plasma velocity, B is the magnetic 
field, P is the gas pressure, g is the gravitational acceleration of a star with mass M* and 
radius i?*, / is the identity matrix, and e is the total energy density, given by 


.^2 P 

Y ^ 7 - 1 ^ Y' 


(5) 


Here, 7 is the poly tropic index such that P (x p'^. 

BATS-R-US has been ehosen for this work as it has been used extensively [|371 |M1 [3^ . 
and reliably models the wind behaviour of the Sun. Additionally, the block refinement within 
the eode allows the resolution to be inereased in speeifie areas of interest, minimising the 
load on computational resources by not having to inerease resolution everywhere. It is a 
modular eode, designed to be used in different ways and with different physieal parameters 
than the original solar parameters. In essenee the code modifications involve the replaeement 
of physical constants such as solar mass and radius sun by stellar parameters that ean be 
varied as needed for eaeh new stellar target. For more details see [|40l HTI . 


Preliminary Results and Discussion 

Some preliminary results from our implementation of the BATS-R-US eode for the star r 
Boo are detailed in Figure The plot shows a shoe in the equatorial plane of the star out 
to 20 stellar radii, with the eolour map indieating the magnitude of the wind velocity in 
km s“^. For context, eontours of radial magnetie field from -0.2 to - 1 - 0.2 Gauss are overlaid 
in blaek (negative field) and white (positive field). The complex magnetie field strueture 
seen in the lobes of positive and negative field eontours produces abrupt differenees in wind 
velocities around the star. Further work is needed to eonfirm that the apparent sharp changes 
in wind veloeity are completely physical, and not an artefact of the limited phase eoverage 
of the magnetie mapping, or the limited resolution of the wind modelling used to produee 
these preliminary results. However, these ‘streams’ of higher velocity wind at the boundary 
between the areas of positive and negative magnetic field are particularly interesting, and 
eould suggest the funnelling of plasma through these boundary regions. 
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Fig. 2: A plot of the wind velocity of the star r Boo from observations taken in July 2008 by 
Fares et al. [l25ll . in the equatorial plane of the star out to 20 stellar radii. Overlaid are contours 
of radial magnetic field ranging from -0.2 to -1-0.2 Gauss, with white indicating positive field 
and black indicating negative field. Showing these contours highlights the interesting ‘streams’ 
of higher wind velocity at the boundaries between positive and negative magnetic field. 


A future goal is to understand the temporal behaviour of the star r Boo, as well as 
understanding the impact its planets will have on the behaviour of the wind. We will explore 
further targets, including others that have been mapped with ZDI, which will be used to 
investigate the time evolution of winds, and some that have not previously been explored 
with ZDI. 

As well as producing wind models for the target stars, we will modify the BATS-R-US code 
to simulate the interaction of the stellar wind with a planet for each target. Beyond allowing 
us to better characterise these particular planetary systems, these simulations will also help us 
to better understand the formation and evolution of planetary atmospheres beyond our Solar 
system. 

It is through that work that we aim to understand the impact of the winds on the planets 
that orbit these stars, and similarly how the presence of a planet, which may or may not have 
a magnetic field of its own, impacts on the behaviour of the wind. 
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